Effects of Line-tying on Magnetohydrodynamic Instabilities and 

Current Sheet Formation 

Yi-Min Huang/'2'3:|] A. Bhattacharjee/'2'3.|t| and Ellen G. Zweibel^'^,! 

^ Space Science Center, University of New Hampshire, Durham, NH 03824 
^Center for Integrated Computation and Analysis of Reconnection and Turbulence 
^Center for Magnetic Self- Organization in Laboratory and Astrophysical Plasmas 
^Department of Physics and Department of Astronomy, 
University of Wisconsin, Madison, Wisconsin 53706 

Abstract 

An overview of some recent progress on magnetohydrodynamic stability and current sheet for- 
mation in a line-tied system is given. Key results on the linear stability of the ideal internal kink 
mode and resistive tearing mode are summarized. For nonlinear problems, a counterexample to 
the recent demonstration of current sheet formation by Low et al. [B. C. Low and A. M. Janse, 
Astrophys. J. 696, 821 (2009)] is presented, and the governing equations for quasi-static evolu- 
tion of a boundary driven, line-tied magnetic field are derived. Some open questions and possible 
strategies to resolve them are discussed. 



* Electronic address: |yimin.huang@unh.edu| 
Electronic address: | amit ava . bhatt achar j ee@unh.edu| 
•l- Electronic address: Iz weibel@astro . wise . edul 



1 



I. INTRODUCTION 



In plasma physics, line-tying refers to the presence of conducting walls enclosing the 
magnetized plasma of interest with a non-vanishing component of the magnetic field normal 
to the wall. Line-tying is usually employed as an idealization of the boundary condition in 
some astrophysical plasmas, where the plasma density varies several orders of magnitude 
over a short distance. For example, the magnetic field lines in the coronae of stars and 
accretion disks are rooted in the dense, highly conducting gas below. In the limit of infinite 
density contrast, the dense gas may be treated as a conducting wall at which the field lines 
are tied. 

Line-tying imposes a strong constraint on the plasma motion, therefore in general is 
stabilizing. This is particularly true in the case of ideal magnetohydrodynamics (MHD). 
In ideal MHD, magnetic field lines are frozen-in to the plasma flow. If the footpoints of 
magnetic field lines are not moving, i.e., the wall is rigid, then the footpoint mapping 
following the field lines from one end to another is conserved. In resistive MHD, magnetic 
field lines are allowed to slip though the plasma on a resistive time scale. In this sense, the 
line-tying becomes imperfect, and the stabilizing effect becomes less stringent. It has long 
been speculated that solar coronal loops can remain stable for longer than typical MHD 
instability time scale due to line- tying stabilization, jl, 2] 

Line-tying may also be used in a broader context with a flow imposed on the boundary, 
to model the convection on the solar surface or the rotation of accretion disks. The imposed 
flow shuffles the footpoints of magnetic field lines, twists up and entangles them, thereby 
converting the kinetic energy of the flow into the magnetic energy in the field. Release 
of the stored magnetic energy powers some of the most splendid plasma phenomena such 
as solar flares and coronal mass ejections (CMEs). Dissipation of the magnetic energy on 
smaller scales is thought to be the energy source of the coronal heating. In a seminal 
paper, Parker argued that when the magnetic field line topology becomes complicated due 
to the footpoint shuffling, tangential discontinuities (current sheets) inevitably will form.j^ 
In the presence of dissipative mechanisms such as resistivity, tangential discontinuities will 
be smoothed out and the magnetic energy will be turned into heat. Parker coined the term 
"topological dissipation" to describe this process. Following Parker's original suggestion, 
many authors 4t9| have shown by direct numerical simulation that random shuffling of 



2 



the coronal magnetic field lines by photospheric motions progressively increases the current 
density, resulting in sporadic energy release. 

The key ingredient of Parker's scenario of coronal heating is the ubiquitous presence 
of current sheets in a magnetized plasma. To illustrate his point, Parker considered a 
uniform magnetic field Bqz line-tied to conducting plates ai z = ±L/2, which represent 
the photosphere. Motions in the photosphere displace the field-line footpoints and entangle 
the field lines. Suppose we freeze the footpoints and let the system relax to the minimum 
energy state under the condition of preserving all the ideal topological constraints. Parker 
argues that in general the minimum energy state must contain tangential discontinuities. 
Since its first appearance in Ref. j^, this simple setting has become the "standard model" 
of the Parker problem. The property that magnetostatic equilibria tend to form tangential 
discontinuities is stated as a "Magnetostatic Theorem" in Parker's book. 10| Although Parker 



gives compelling physical arguments in support of the Theorem, no rigorous proof is given. 

Parker's claim has stimulated considerable debate that continues to this day ll-[l7|, with- 
out apparent consensus. The point of contention is whether or not a magnetic field with 
neither null points, nor separatrices, nor closed field lines can develop tangential discon- 
tinuities. It is well established that null points, separatrices, or closed field lines are the 
preferential sites for current sheet formation when a system undergoes some instabilities or 



is deformed via boundary conditions. 



18|-|24|| In the Parker problem, however, none of these 



structures are present. Van Ballegooijen was the first to question Parker's claim. He demon- 
strated by an analytical argument that the equilibrium will always be continuous, unless the 



footpoint mapping is discontinuous. 
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25l | Subsequently, several authors have found smooth 
equilibria using either linear theories for small footpoint displacements or Lagrangian nu- 



11, 



26| However, as pointed out 



merical relaxation schemes for nonlinear problems. |12l . 

by Ng and Bhattacharjee, 15|] van Ballegooijen's argument is valid only for simple current 
sheets; it fails if a current sheet has multiple branches joined together (see Figure 2 of Ref. 



Ng and Bhattacharjee also proved a theorem which asserts that there is at most one 
smooth equilibrium for any given smooth footpoint mapping. It follows that if a smooth 
equilibrium becomes unstable, there is no other smooth equilibrium the system can relax to; 
therefore the system must relax to an equilibrium with tangential discontinuities. The the- 
orem of Ng and Bhattacharjee is within the framework of reduced magnetohydrodynamics 
(RMHD), which assumes the existence of a strong guide field. [27|, l28| It also assumes double- 
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periodicity in the plane perpendicular to the guide field. It is not known at the present time 
whether the result can be generalized to full MHD or more general boundary conditions, 



although Aly 29| appears to have taken an important step towards that goal. Lagrangian re- 



laxation studies indeed suggest current sheet may form under certain circumstances. |16l. |17| 
although it is rather difficult to draw decisive conclusions from numerical studies due to 
limits of spatial resolution. 

This paper gives an overview of our recent progress on MHD stability and current sheet 
formation in a line-tied system. Previously published results are summarized in Sec. [IT] and 
Sec. Illlt with new perspectives added. Sec. |TT] summarizes the results on linear stability 
regarding line-tying effects on the ideal internal kink mode and resistive tearing mode. 
Parker's problem on current sheet formation is discussed in Sec. IIIII The discussion is 
divided into two parts. The first part discusses why the nonlinear evolution of ideal kink 
mode is a relevant setting to test Parker's hypothesis. The second part is a discussion 



about the recent demonstration of current sheet formation by Low et al. 30|-|33|] and our 



objection to it. 3J] Some outstanding open questions will be discussed along the way. Sec. 
IIVI contains the new results, where the governing equations for quasi-static evolution of a 
boundary driven, line-tied magnetic field are derived, and a possible strategy of applying 
these equations to fully settle an open question regarding the Parker problem as well as the 
demonstration of Low et al. is outlined. Finally, some future perspectives are discussed in 
Sec. IVl 

II. LINE-TYING EFFECTS ON IDEAL AND RESISTIVE MHD INSTABILITIES 
A. Ideal Internal Kink Mode 

A line-tied screw pinch is often employed as a model of solar coronal loops. Assuming 
that the plasma pressure is negligible, in cylindrical coordinates (r, 0, 2;), the equilibrium 
magnetic field is 

B = 5^(r)0 + 5,(r)z, (1) 

which satisfies the force balance equation 

d 52 



0. (2) 
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In a periodic pinch of length L, we may decompose the plasma displacement $, into Fourier 
basis in both z and (p directions. Therefore we consider an eigenmode of the form ^ = 
^(^^^g-yt+i(kz+m(f)) ^ where k = 2mi/L, 7 is the growth rate, and both n and m are integers. 
The m = ±1 internal kink mode is unstable if the safety factor q = 2TxrBz/ LB^ drops below 
unity within the plasma and increases with radius. Unstable modes have helical symmetry, 
depending on z and in the combination kz ± 0, and possess a so-called resonant surface 
r = Tg, on which k - B = kEz^B^/rg = 0. Due to the constraint of periodicity, the resonant 
surface must also be a rational surface on which magnetic field lines are closed. The driver 
of the instability lies within r^, and the radial component of the plasma displacement nearly 
vanishes outside it. It has been shown that the thin transition layer near predicted from 
linear theory corresponds to an infinite current sheet in finite amplitude theory, at least 
within the framework of reduced, ideal MHD. 19|, |35|] In a line-tied pinch, the boundary 
condition constrains the plasma displacement $, to zero at z = ±L/2. Therefore a Fourier 
decomposition can no longer be done in z, and the eigenvalue problem is inherently two- 
dimensional (2D). At a fundamental level, there is a deeper issue. In a line-tied screw pinch, 
all magnetic field lines are topologically equal; they all go from one end plate to another 
(assuming ^ everywhere). Therefore the notions of rational surfaces and resonant 
surfaces are no longer valid. A question is: is there still a steep gradient in the eigenmode, 
as in the periodic counterpart? And if the answer is yes, where will the steep gradient be? 



In Ref. 361] a detailed study of the line-tied internal kink mode has been carried out. 
It was found that the fastest growing internal kink mode in a line-tied system possesses 
a steep internal layer. The internal layer is located at the resonant surface of the fastest 
growing periodic mode in an infinitely long system (i.e., k is a continuum). However, the 
internal layer in a line-tied system is smoother than its periodic counterpart. Only in the 
limit L -> 00 does the former approaches the latter. This can be seen in Figure [H which 
shows the radial displacement and parallel current Jii of eigenmodes for different L along 



0, as compared to its periodic coun 



;erparts. Similar results have been 



the midplane z 

found in Refs. 37|, l38| for different equilibria and Ref. (39| including finite pressure effects. 

There is a critical length for the instability. When L < L^. the internal kink is 
stabilized. The critical length Lc is approximately given by 

U^2\Va.I/^o, (3) 
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where 70 and are the growth rate and the resonant surface of the fastest growing periodic 
mode, respectively, and Vaz is the z component of the Alfven speed. This is roughly con- 
sistent with the heuristic stability criterion that the unstable mode is stabilized when the 
Alfven transit time along the loop is less than the inverse of the growth rate in a periodic 
system. Once the critical length is known, the growth rate 7 for a system of length L 
can be determined by solving the equation 

^ = Ip^ . (4) 

L tanh ^(7/70) 

In the limit L — j- 00, the growth rate 7 — ?■ 70. The growth rate 70 in a periodic pinch scales 

as 

7o~e'|V:4.|/a, (5) 

where a is the characteristic radius of the pinch and e ~ B^p/Bz is the tilt angle, which 
typically is a smaller parameter. The critical length Lc ^ a/e^. In solar applications, often 
the coronal loop length and the axial field Bz are given, and the azimuthal component B^ 
is generated by applying some footpoint twisting. Then we can ask the following question: 
what is the critical twisting angle (pc for the loop to become unstable? From ~ a/e'^ we 
have the critical tilt angle ~ (a/LY^'^ when L is given. That corresponds to a critical 
twisting angle 



a \a / 



(6) 

The precise value of the critical tilt angle ec depends on the profile. For the profile i?^ = 



const X (r/a) exp(— (r/a)^), which was considered in Ref. 36|, the critical tilt angle at r = a 
is approximately given by 

e,^l.lx(^-j . (7) 

The critical tilt angle depends on the aspect ratio (a/L) rather weakly. For a coronal loop 
with a/L = 1/100, the critical tilt angle ec — 0.24, which is approximately 14°. For a shorter 
loop with a/L = 1/10, the critical tilt angle is about 29°; whereas for an extremely long loop 
with a/L = 1/1000, the critical tilt angle is about 6.3°. These numbers are in reasonable 



agreement with estimates obtained from balancing the Poynting flux from the coronal 



and the observed coronal loss rate, which give a tilt angle ranging from 10° to 20°. 40 



base 



4l| 
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B. Resistive Tearing Mode 



When a small but finite resistivity r] is included in the theory, line-tying becomes im- 



perfect. In Ref. 



421 ]. Delzanno et al. studied the resistive effects on line-tied internal kink 



modes. They found that when the system length L is greater than the ideal critical length 
jjideai ^ the mode growth rate is similar to the ideal one. However, resistivity can make a 
system with L < L^f'^"-'' unstable. Interestingly, the critical length of a resistive system LJ!*^* 
is found to be independent of r]. 



In Ref. 



43| , we studied the effects of line-tying on the resistive tearing instability in slab 



geometry. The equilibrium has a strong guide field along z, and By changes direction 
across a current layer around x = 0. Two conducting plates a,t z = ±L/2 provide the line- 
tied boundary condition. This configuration is ideally stable, but can be unstable to the 
tearing mode when rj > 0. The instability is studied within the framework of the RMHD 



equations. 
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441 ] After proper normalization, the magnetic field and the velocity can be 



expressed in terms of the flux function ip and the stream function as B = z + x z 
and u = Vj,^ x z. The governing equations are 



dtn + [<p,n] = d,j + [^,j], (8) 

Here V_l = + ydy is the perpendicular gradient; Q = —V^cp , J = —V^ip are the 
vorticity and the electric current; and [ip, </>] = dyipdx4> — dxipdy(f) is the Poisson bracket. The 
standard Harris sheet By = const x tanhx, is assumed, but more general profiles can be 
easily implemented. 

In a periodic system, we may consider linear perturbations of the form (p = 
0(x)e^*"'"*'^f^"'"*'^'^^ and ip = ■?^(x)e'^*"'"*'^f^"'"*'^^^. The system is unstable to the tearing mode 
with a growth rate 7 ~ 77^/^ in the asymptotic limit 77 — 0. With line-tying, Fourier de- 
composition along z is not possible and a two-dimensional eigenvalue problem has to be 
solved. In this case, the system is unstable only when L is greater than a critical length Lc, 
which is found to be independent of rj. However, near marginal stability, the linear unstable 
mode grows slowly, and plasma inertia can be neglected everywhere. With this "force-free" 
approximation, the growth rate 7 is proportional to rj. When L is sufficiently long, plasma 
inertia can no longer be neglected within the resistive layer. In this regime, the line-tied 
growth rate approaches the periodic one, and the 7 ~ 77^/^ scaling is recovered. The transi- 
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tion from the force-free regime to the inertia dominated regime occurs at a transition length 
Lt which scales as Lt ~ ?7~^/^. Notice that without line- tying, the force- free approximation 
only applies in the outer region; in the resistive layer, plasma inertia (or viscosity) has to be 
included, otherwise the equations become singular at the resonant surface where k ■ B = 0. 
In a line-tied system, on the other hand, it is possible that the force-free approximation is 
applicable over the whole domain, provided that the mode grows slowly. This is because 
that there is no resonant surface in a line-tied system, and the equations are not singular 
even without plasma inertia. The resistive layer persists even with line-tying. However, line- 
tying effects make the resistive layer smoother; when the system length L becomes longer, 
the resistive layer becomes steeper. This general trend, as shown in Figure [21 resembles 
that in the internal layer of the ideal internal kink mode. Similar results have been found 
in cylindrical geometry by Delzanno et al. 45| 



III. CURRENT SHEET FORMATION AND THE PARKER PROBLEM 

A. Nonlinear Saturation of the Ideal Kink in Line-tied Geometry 

A question that remains open is the nonlinear saturation of ideal internal kink mode in 
a line-tied screw pinch. After the onset of the instability, if we let the magnetic field relax 
to a minimum energy state but preserve all the ideal topological constraints, will it develop 
a tangential discontinuity? This is an excellent setting for testing Parker's conjecture, as 
the initial condition contains neither null points nor separatrices nor closed field lines. In 
a periodic pinch the answer to the question is yes, as shown by Rosenbluth et al. that a 
helical ribbon-like current sheet will form at the resonant surface, which is also a rational 
MU a .i„e-t.d .yste,„. the ..we. . .ot so e.ea. The steep g.adie„t the 



surface. 



line-tied eigenfunction (Figure d]) suggests that a thin current filament will form. Indeed 
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a helical thin current filament has been observed in several nonlinear simulations. 
However, whether it is a current singularity or not is very difficult to settle by a numerical 
calculation. We have seen that line-tying in general has the effect of smoothing out the 
internal layer of the linear eigenmode; it may have a similar effect to the nonlinear saturation 
and turn a current singularity into a thin current layer with a finite width. If this turns 
out to be the case, we may expect the current layer to become thinner as the system length 
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becomes longer, assuming that the equihbrium profile remains unchanged. 

The theorem of Ng and Bhattacharjee, as discussed in the Introduction, asserts that 
if a smooth equilibrium becomes unstable, the system must relax to an equilibrium with 
tangential discontinuity. According to the theorem, it appears that the ideal internal kink 
mode will develop a current sheet. However, the theorem does not apply to the present case 
for the following reasons. First, the theorem is proved within the framework of RMHD, in 
which the ideal internal kink mode is marginally stable. To capture the internal kink mode 
one has to carry out the expansion to higher orders. This also raises the question whether 
the theorem can be generalized to full MHD or not, as some unstable modes are clearly 
missing in RMHD. Second, the proof assumes double periodicity in the perpendicular plane, 
and a screw pinch does not belong to this category. 



B. A Counterexample to the Demonstration of Current Sheet Formation by Low 
et al. 



Recently, in a series of papers 



30- 



Low and coauthors have tried to demonstrate 
unambiguously the formation of current singularities using what they call "topologically 
untwisted fields", which, in the context of force-free fields, is synonymous with potential 
fields. That is, they limit themselves to a special subset of force-free fields which contains 
no electric current. In this case, the magnetic field can be expressed as B = Vx with 
some potential Xi where x satisfies V^x = because B is divergenceless. In a simply 
connected, compact domain, x is uniquely determined up to an additive constant if the 
normal derivative of x is specified on the boundary. Therefore, a potential field is uniquely 
determined by prescribing the normal component of B (= n- Vx, where n is the unit normal 



vector) on the boundary. In Ref. (32|, the authors consider a potential field in a cylinder 
of finite length as an initial condition; the normal component of the field is nonzero only 
at the top and the bottom of the cylinder. The cylinder is then compressed to a shorter 
length. Because the whole process is governed by the ideal induction equation, the normal 
component of the field remains the same. By making a key assumption that the field will 
remain potential during the process, they circumvent the complicated problem of solving 
for the various stages of quasi-static evolution, and simply calculate the magnetic field in 
the final state as the potential field which satisfies the boundary conditions. They then 



9 



numerically calculate the magnetic field line mapping from one end to the other, for both 
the initial and the final states. They find that for a three-dimensional (3D) field, the field 
line mapping, and hence, the topology, is changed. Because an ideal evolution preserves the 
magnetic topology, the authors claim that the field must evolve into a state which contains 
tangential discontinuities; furthermore, singularities may form densely due to the ubiquitous 
change of the field line mapping. 

The assumption that the field will remain potential during the compression, however, 
is not proven in their work. The reasoning behind the assumption could be understood 
as follows. Consider a magnetic field line (the base field line) and a neighboring field line 
which is infinitesimally separated from it. As we follow along the field lines, in general the 
neighboring field line winds around the base field line at a rate d(j)/ds 0, where s is the arc 
length along the field line and (j) is the azimuthal angle of the displacement vector between 
the two field lines relative to a comoving reference frame. Here the comoving reference frame 
is specified by three mutually perpendicular vectors at each point along the field line. One 
of the vectors is tangent to the field line, while the other two follow parallel transport along 
the field line. (See Appendix |A] for the detail of how parallel transport is defined. One also 
should not confuse the azimuthal angle (p here with the angle of rotational transform in a 
toroidal device, where the base field line is a closed loop. For the rotational transform, the 
well-known Frenet-Serret frame is a good reference frame. The Frenet-Serret frame does not 
follow parallel transport in case the base field line has a nonvanishing torsion. See Appendix 
|A]for further discussion.) In general, different neighboring field lines have different winding 
rates. However, for a force-free field satisfying 



with B ■ VA = 0, it can be shown that the averaged local winding rate over all neighboring 
field lines is related to A as (see Appendix |A] for a derivation) 



This simple relation shows that parallel current is related to the winding between neighboring 
field lines. For a potential field, which satisfies A = everywhere, the average local winding 
rate {d(f)/ds) is equal to zero. This justifies the use of the term "untwisted fields" for 
potential fields. Intuitively, it would appear to require a nonzero vorticity on the boundary, 



V X B = AB 



(9) 




(10) 
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which is absent in a simple expansion or compression, to twist up the field. Therefore, 
it is plausible that a potential field will remain potential during a simple expansion or 



compression. However, in Ref. [3J], an explicit counterexample is constructed in a two- 
dimensional (2D) slab geometry. The counterexample demonstrates that a potential field 
can evolve to a smooth, non-potential field under simple expansion or compression. The 
construction is summarized as follows. 

Consider a 2D configuration in Cartesian geometry, with x the direction of symmetry. A 
general magnetic field may be expressed as 

B = B^x + X X V?/', (11) 

where = Bx{y, z) and ip = ip{y, z). For a force- free field, = Bx{ip) and ip must satisfy 
the Grad-Shafranov equation 

V^i. ^ -B,^, (12) 

where A in Eq. ((21) is equal to —dBx/dil) in this representation. If the field is potential, then 
Bx = const and V'^ip = 0. 

Let us consider a force-free field bounded by two conducting surfaces at z = and z = L, 
with all field lines connecting one end plate to the other. The footpoint mapping from one 
end to the other is characterized by tp on the boundaries, as well as the axial displacement 
along X when following a field line from the bottom to the top: 

r-Z = L 

-0 



h{^) = B^iij) { I T^ds] . (13) 



Here ds is the line element on the y — z plane, and the subscript ijj indicates that the 
integration is done along a constant ip contour. 

Let the initial condition be (-B^O) V^o) with a system length Lq. Suppose the system has 
undergone a simple expansion or compression such that the system length changes to some 
L Lq. Since ideal evolution preserves the footpoint mapping, the fiux function ip on the 
boundaries and the axial displacement h{ip) must remain unchanged. To determine the final 
force-free equilibrium, we must solve the two coupled equations 

V^. ^ -ftf , (14) 

and 
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subject to the conditions ip{y,0) = ipo{y,0), ip{y,L) = ipQ^y , Lq) , where the axial displace- 
ment h{il)) is determined from the initial condition. The set of coupled equations fll4l) and 
( IT5|) is called a generalized differential equation, which in general requires numerical solu- 



tions. 



52| 



In Ref. j3^, we consider the initial condition t/^o = y + e sin(y) cosh(2 — 1/2), B^q = 1, and 
Lq = 1, where e -C 1 is a small parameter. Eqns. f[T^ and f[TB]l are solved for arbitrary L both 
analytically with a perturbation theory and numerically. The results from both approaches 
are consistent, and the final state is found to be smooth, but contains non- vanishing current, 
and is therefore non-potential. 

It should be pointed out that our configuration is not within a compact domain. Further- 
more, the example we worked out is periodic along the y direction, even though that is not 
required by the formalism. Topologically, the geometry is equivalent to the space between 
two concentric, infinitely long cylinders, therefore not simply connected. 53|] Consequently, 
specifying the normal component of the magnetic field on the boundary does not uniquely 
determine the potential field. One can, for example, add a constant By or B^ to the field 
without changing the normal component of B. In this regard, our configuration does not 
belong to the same class of configurations considered by Low et al. Nevertheless, our exam- 
ple shows that the demonstration by Low et al. is not complete if their key assumption is 
not proven. Another criticism of the demonstration of Low et al. has been raised recently 



by Aly and Amari from a different point of view.fSJ] To fully settle the issue in the future, 
a more general formalism for the quasi-static evolution of a boundary driven magnetic field 
is needed. This is the topic of the next section. 

IV. QUASI-STATIC EVOLUTION OF A BOUNDARY DRIVEN LINE- TIED 
MAGNETIC FIELD 

The governing equation for quasi-static evolution driven by boundary motions was derived 



by van Ballegooijen (VB) Hj within the framework of RMHD. This is under the assumption 
that the characteristic time scales of the boundary motions are much longer compared to 
the Alfven transit times, therefore the system remains force-free for all time to a good 
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approximation. Here we generalize the VB equation to full MHD without assuming any 
asymptotic ordering. 

We assume that the system satisfies the force-free condition for all time. The magnetic 
field evolves according to the ideal induction equation 

= -V X E = V X (u X B) , (16) 

where u is the plasma fiow and E the electric field. Now suppose the magnetic field B and 
the boundary velocity are given at t = to- The objective is to find a plasma fiow u which 
is consistent with the prescribed boundary fiow u^, and which, when evolving the magnetic 
field according to the ideal induction equation (ITB]) . carries the system to a neighboring 
force- free equilibrium at t = to + dt. Following VB, first we take the time derivative of the 
force-free condition ([9]). This yields 

VxdtB = BdtX + XdtB. (17) 

Using the induction equation (fT6ll and B ■ VA = 0, Eq. ( |T71) may be rewritten as 

— B = (9tA + u ■ VA) B = V X (AE - V X E) 

dt 

= -V X (Au X B- V X (u X B)). (18) 

Taking the divergence of Eq. (ITSl) yields 

B.v|^0, (19) 

which is consistent with the condition that B ■ VA = is satisfied for all time. Equation 
( TTSll is the equivalent of the VB equation in the framework of full MHD. It has to be solved 
for u subject to the boundary condition prescribed by the boundary fiow u''. In other 
words, we have to find a fiow u such that V x (Au x B — V x (u x B)) is proportional to 
B everywhere inside the domain. Clearly only the component of u perpendicular to B, u^, 
can be determined from this approach. The parallel component u\\ is arbitrary, except that 
it has to be consistent with the boundary condition. Equation (1181) may be solved for very 
general boundary conditions, including convection (u^ ■ n = 0), expansion or compression 
(u* X n = 0), or a combination of both. If the trajectory of each boundary point is prescribed, 
solving Eq. (1181) at each instance and integrating Eq. (1161) at the same time will carry the 
magnetic field through a series of force-free equilibria. Notice that Eq. flTSl) is linear in 
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u at each instant of time. That is, if {dXi/dt;ui) is a solution for the boundary flow u\ 
and {dX2/dt;u2) a solution for then {cidXi/dt + C2dX2 / dt , CiUi + C2U2) is a solution 
= CiuJ + C2U2, where Ci and C2 are arbitrary constants. 

The time evolution of A is of particular interest here. It allows us to address the issue 
regarding Low's assumption. Starting from a potential field, i.e., A = everywhere, Eq. 
( ITS]) becomes 

= V X V X (u X B) . (20) 

The field will remain potential if dX/dt = everywhere. Therefore, it will be sufficient 
to take the initial condition of Low et at, and solve Eq. fl20p . The boundary conditions 
can be chosen as u^^^ = z, u^^^^^^ = 0, and u^^^ varies continuously from the bottom to 
the top and matches the velocity at the bottom and the top. If the solution turns out to 
have dX/dt 7^ 0, then the assumption of Low et al. is disproved. On the other hand, if 
Eq. fl2Ul) does not have a smooth solution, then current singularities will form immediately. 
A third possibility is that a smooth solution exists, and satisfies dX/dt = 0. In this case 
we have to follow the time evolution through a series of potential fields until either one of 
the previous two possibilities is realized. (The field cannot remain a smooth potential field 



forever, as this contradicts the findings of Ref. |33|.) If Low et al. are indeed correct, this 
approach may answer the questions regarding when, and where current singularities will 
form. It is not difficult to see that in general dX/dt will not be zero, as follows. The solution 
to Eq. fl2Dl) determines the components of u perpendicular to B; that is, there are two 
degrees of freedom at each point. And Eq. fl2Up demands that V x V x (u x B) be parallel 
to B, which again corresponds to two independent constraints at each point. Therefore 
we have the same number of unknowns and equations. If the magnetic field is to remain 
potential, then the condition (V x V x (u x B))|| =0 should be satisfied as well. This 
would correspond to a third independent constraint at each point, and the system is then 
overdetermined. Of course, the boundary motion considered by Low et al. is a very special 
one. Therefore we cannot preclude the possibility that in their case V x V x (u x B) || B 
implies V x V x (u x B) = 0, when the solution is smooth. If this can be proven to be the 
case, then the assumption of Low et al. is proven. 

The approach outlined here has been applied to the example described in Sec. IIII B[ 
And indeed the solution of Eq. (1201) has dX/dt 7^ 0, consistent with our findings from the 
approach mentioned in Sec. IIII Bl that the field will evolve into a non-potential field. A 
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similar approach has been taken by Scheper and Hassam in their demonstration of current 



2J] The initial 



sheet formation in 2D slab geometry, for an initial field with separatrices. 
condition of Low et al. is fully 3D, therefore solving Eq. (1201) is more difficult, requiring 
numerical treatment. 



V. FUTURE PERSPECTIVES 

In this paper, we have reviewed some of the recent progress regarding the effects of line- 
tying on MHD instabilities and current sheet formation and presented some new results. For 
the problems of linear stability, it is shown that line-tying in general has a stabilizing effect. 
Furthermore, because of the absence of a resonant surface in a line-tied system, the internal 
layer of the eigenfunction, although it persists, is smoother than its periodic counterpart. 
For the more important nonlinear problems, there remain two outstanding questions: the 
nonlinear saturation of the line-tied ideal internal kink mode and the complete resolution 
of the issue regarding the assumption made by Low et al. in their demonstration of current 
singularities. These open questions pose significant challenges for the future, involving 
analysis as well as numerical computation. The nonlinear saturation of the internal kink 
mode may be attacked by some numerical force-free field solvers that preserve the field line 



topology, such as Lagrangian relaxation schemes. llTl . |55| If the final equilibrium has a 
thin, but non-singular current layer, it may be fully resolved with convergence, provided 
that the resolution is sufficient. For the second question, a possible approach is outlined 
in Sec. llVt by solving the derived governing equations for a quasi-static, boundary driven 
magnetic field. Another possible approach is solving the final equilibrium of the compressed 
system by a force-free field solver, mentioned above. In any case, we are facing a problem 
where the solution may become singular, therefore a reliable mean to ascertain singularities 
numerically may be needed. 

The formalism derived in Sec. IIVI also opens up the possibility of generalizing the proof 
of Ng and Bhattacharjee, which relies heavily on the VB equation, to full MHD. This will 
be a topic of future investigation. 
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Appendix A: Comoving Frame of Reference, Parallel Transport, and a Derivation 
ofEq. ([To]) 

Let X(s) be a magnetic field line (the "base" field line) parameterized by the arc length 
s. A comoving reference frame along the magnetic field line is specified by three mutually 
perpendicular unit vectors {ei(s), 62(5), 63(3)} at each point along the field line. Here we 
further require 63 to be the tangent vector. That is, = d^/ds = b = B/i?. In general, 
the basis vectors {Gi}-^^ 23 follow 

dei 

— ^ = cj X ei, (Al) 
ds 

where u: can be regarded as the angular "velocity", if we treat the arc length s as the time 
variable. We adopt the following definition of parallel transport: a comoving frame is said 
to be parallelly transported along a field line if the angular velocity is perpendicular to the 
magnetic field. That is, the angular velocity satisfies u ■ = 0. Under this condition, the 
comoving frame is not rotating with respect to the direction of the magnetic field, therefore 
is a good reference frame for measuring the local winding rate of neighboring field lines. 
From de^/ds = u; x 63 and u; ■ 63 = 0, the angular velocity is uniquely determined as 

. = e3 X -J, (A2) 

If a set of basis vectors is given at one point, then the comoving frame following parallel 
transport at each point can be obtained by integrating Eq. (lAlj) with the angular velocity 
given by Eq. I \A2\} . Specifically, 

dei ( deA 
de2 ( c/egN 

d7 = -r^- drr" 
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Having set up the comoving frame, we can now define a local coordinate system (oi, 02, s) 
in a neighborhood of the field line: 



x(ai, a2, s) = X(s) + aiei + 0262. 



(A5) 



Now suppose there is a neighboring field line with a trajectory defined by X'(s) = X(s) + 
5x(s) = X(s) + ai(s)ei + 02(5)62, where ai and 02 are infinitesimal. The differential arc 
length of the neighboring field line X'(s) is 



ds' 



IdX.' c/X' 



ds = ( 1 — 5x ■ — - I c/s, 
V ds 



ds ds 

where Eqs. (1A3I) and (1A4I) are apphed. The magnetic field at X'(s) is 

B'(s) = S(s)e3 + (5x ■ VB, 

which has a magnitude 

= + 63 ■ ((5x ■ VB) . 
Using Eqs. (|A6l), dAT]), and (IA8]) in the field hue equation d'X.'/ds' 



(A6) 



governing equation for (5x: 



(5x 



((5x ■ vb; 



(AT) 
(A8) 

B'/B' yields the 
(A9) 



\ds "J ^ B{s) 

where the subscript _L indicates a projection onto the plane spanned by ei and 62. Equation 
(1A9I) can be written explicitly as 



d_ 

ds 





= M{s) 


ai 


02 




0-2 



(AlO) 



where 



and 



M{s) 



Bis] 



[VB)n (VB)2i 

[VB)i2 (VB)22 



[VB)i,- = (e,- VB)-e,-. 



(All) 



(A12) 



The winding rate of the neighboring field line becomes apparent if we express the deviation 
5x in terms of the local polar coordinates (r, 0), which are related to (ai, 02) by ai = r cos 
and a2 = r sin 0. It follows from Eq. flAlOp that 
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Idr 1 , (VB)i2 + (VB)2i . , . 2 



r (is -B(s) 

and 



(VB)ii cos' + ^ '-^^—^ — sin 20 + (VB)22 sin' , (A13) 



[VB)22-(VB)n . , 2^ ,w-oN ■ 2 



■ sin 20 + (VB)i2 cos' - (VB)2i sin' . (A14) 



B{s) \ 2 

Equation ( ]A13[) determines whether the two field fines are diverging from or converging 
towards each other, and Eq. (lAMD gives the winding rate d(j)/ds of the neighboring field 
line with respect to the base field line. Clearly, d(j)/ds depends on 0. Therefore, in general 
different neighboring field lines have different winding rates. Averaging Eqs. (1A13P and 
f lA14p over yields 



and 



((VB)i2 - (VB)2i) = (A16) 



\dsl 2B{s)'' '^^ ' IB 
Here we have made use of the condition V ■ B = in the last step of Eq. (lAlSp . Equation 
(lAlSP corresponds to the conservation of magnetic flux, d{(jB)/ds = 0, where a is the cross 
section of a thin flux tube. And Eq. flA16P relates the average winding rate to the parallel 



current. For a force-free fleld satisfying V x B = AB, Eq. (]A16P can be further simplified 

as 

and we have thus completed the derivation of Eq. fITU]) . Heuristically, Eq. flAlSp may be 
understood by consider a circular cross section of a radius r. The conservation of magnetic 
fiux, d{aB)/ds = 0, implies a~^da/ds = —B^^dB/ds. Using this equation with the relation 

.da 1 / dr , / 2 dr\ , , . 

a-'— = —6 27rr-rf0 = ( ) , A18 
as Trr^ / as \r as 



Eq. (lAlSP is obtained. Likewise, Eq. (]A16P may be understood by considering a loop 



integral around the circular cross section and using the Stokes' theorem, which yields 

^ rB^dcf) = 27rr {B^) = Trr'eg ■ V x B. (A19) 
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Substituting Eq. flXToD in 

ld(t)\ (B^) 



(A20) 



\ds / rB 

Eq. flA16p is recovered. Our derivation puts Eqs. flAlSp and flA16p on a rigorous footing. 

It is of great interest to consider the well-known Frenet-Serret comoving frame in con- 
junction with the present formulism. The Frenet-Serret frame is defined by the tangent 
vector §3 = (iX/ds, the principal normal vector ei = \de3/ds\~'^ (de^/ds), and the binor- 
mal vector §2 = 63 x ei. These vectors satisfy the Frenet-Serret relations de^/ds = kBi, 
dei/ds = —ne^ + re2, and de.2/ds = — rei, where k, is the curvature and r is the torsion of 



the magnetic field line. 561] The Frenet-Serret relations can be expressed as 

-^Bi =UJ XBi, (A21) 

ds 

with a) = ^§2 + re^. Because u ■ ^ 0, the Frenet-Serret frame does not follow parallel 
transport. The angular velocity of the the Frenet-Serret frame, u), is related to the parallel 
transport angular velocity, a;, by a; = oj + re^. Therefore, the Frenet-Serret frame is rotating 
with respect to a parallelly transported comoving frame with an angular velocity re^. If we 
instead measure the winding rate of a neighboring field line by the angle with respect to 
§1 of the Frenet-Serret frame, then 

f = ^-r. (A22) 
ds ds 

In a toroidal device with its magnetic axis being a closed base field line, we are interested 
in the angle between the initial and final 5x after a full toroidal circuit. The average of this 
angle over many circuits leads to the notion of rotational transform. To measure this angle, 
we need a reference frame that goes back to itself after a full circuit, a requirement that in 
general is not satisfied by a parallelly transported comoving frame. The Frenet-Serret frame, 
on the other hand, only depends on the local geometry of a curve, therefore is well-suited 
for the purpose. As such, the rotation angle after a full circuit can be expressed as a loop 
integral j>d(f) = dcj) — t ds ^ which has contributions from two parts. The first part, ^ dip, 
depends on the parallel current, as discussed above. However, the second part, — <^t ds, only 
depends on the geometry of the magnetic axis. Consequently, a toroidal device can have 
a nonvanishing rotational transform even with a vacuum field, provided that the magnetic 
axis is nonplanar (therefore r 7^ 0). This geometric angle is the basis of the figure-eight 
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and other stellarators with nonplanar axes. Interested readers are referred to Ref. 57|] for 
further discussion. 
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Figure 1: The radial displacement and parallel current Jy of eigenmodes for different L, along 
the midplane z = 0. Also shown for comparison is the periodic case. 
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Figure 2: (Color online) The linear cut of the perturbed stream function (p along the mid-plane 
z = 0, for line-tied cases. Also shown for comparison is the periodic case with kz = 0. 
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